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On the occasion of the 70's birthday of Prof. Adriano Di Giacomo we report on recent numerical 
I ' computations of the Landau gauge gluon and ghost propagators as well as of a non-symmetric MOM- 

scheme ghost-gluon vertex in quenched and full lattice QCD. Special emphasis is paid to the Gribov 
copy problem and to the unquenching effect. The corresponding running coupling a s (q 2 ) is found 
and shown to decrease for q 2 < 0.3 GeV 2 in the infrared limit. No indication for a non-trivial 
infrared fixed point is seen in agreement with findings from truncated systems of Dyson-Schwinger 
equations treated on a four-dimensional torus. 
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1 Introduction 



Three of us (E.-M.I., M.M.-P. and A.S.) had the occasion to meet Adriano Di Giacomo 
many times since the late eighties. Over the years we shared with him the scientific 
interest on topological aspects in connection with the confinement property of Yang- 
Mills theories. For us it was always helpful and exciting to listen to and to discuss with 
him and to follow his ideas and his opinion. We would like to wish him all the best and 
to be continuously active in science in the future. 

In our contribution to this honorary collection of papers we would like to report on 
a field theoretic investigation closely related to the confinement problem which became 
popular over the last decade. It considers the infrared behaviour of the Landau gauge 
gluon and ghost propagators in QCD in connection with the Gribov-Zwanziger horizon 
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condition and the Kugo-Ojima criterion . The propagators also provide a nonpertur- 
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bative determination of the running QCD coupling a s (q ) in a MOM scheme. It is a 
great advantage to be able to compute the relevant Green functions within the framework 
of truncated systems of Dyson-Schwinger (DS) equations as well as within the lattice 
approach and to check the consistency of both approaches. As we shall discuss below, 
for the infinite volume case there is still a mismatch between the infrared limit results, 
what requires further investigations in the near future. 

In the infinite 4- volume limit the DS approach has led to an intertwined infrared power 

behaviour of the gluon and ghost dressing functions ^ 1 1 1 11 

Z D (q 2 ) = q 2 D(q 2 ) oc (, 2 ) 2 \ 

Z G {q 2 ) = q 2 G(q 2 ) cx (q 2 )- X (1) 

with the same A m 0.59 , i.e. a vanishing gluon propagator D(q 2 ) in the infrared occurs 
in intimate connection with a diverging ghost propagator G(q 2 ). Under the condition 
that the ghost-gluon vertex renormalization function Zi(/j, 2 ) is finite and constant (seeEl 



for the perturbative proof ancflSlfor a first SU(2) lattice study) the corresponding MOM 
scheme running coupling is determined by 

a s (q 2 ) = |1 Z D {q 2 ) Z 2 (q 2 ) . (2) 

j— i Isl 
This means that Eq. (UJ) leads to a non-trivial fixed point of a s in the infrared limit? 2 '. As 

we shall show below, our lattice simulations do not agree with this prediction. The strong 

coupling a s (q 2 ) seems to tend to zero in this limit as lattice results for other definitions 

of the running coupling also indicate (se c ) ^^^ for the three-gluon vertex and^l for the 
quark-gluon vertex). 

Possible solutions of the discrepancy have been discussed from different points of 

view in the recent literature. InEEl the DS approach has been studied on a finite 4- 
torus with the same truncated set of equations as for the infinite volume. a s (q 2 ) was 
shown to tend to zero for q 2 — > in one-to-one correspondence with what one finds 
on the lattice. This would indicate very strong finite-size effects and a slow convergence 
to the infinite-volume limit. However, on the lattice we do not find any indication for 
such a strong finite-size effect, except the convergence to the infinite- volume limit would 
be extremely slow. An alternative resolution of the problem has been proposed by Bou- 

caud et alM^.. These authors argued the ghost-gluon vertex in the infrared might contain 
q 2 -dependent contributions which could modify the DS results for the mentioned propa- 
gators a . However, recent detailed DS studies of the ghost-gluon vertex did not provide 

hints for such a modification! 1211SI Thus, at present there seems to be no solution of the 
puzzle. 

For SU{2) extensive lattice investigations can be found inli^l Unfortunately, the 
authors did not reach the interesting infrared region, where the mentioned inconsistencies 
become visible. The same holds for previous SU(3) lattice computations of the ghost and 
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gluon propagators as reported ii i 1 1 I We have been pursuing an analogous study for 
the SU(3) case with special emphasis on the Gribov copy problem^!! Moreover, we 
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have investigated the spectral properties of the Faddeev-Popov operator 1 ^-* 1 . The low- 
lying eigenmodes of the latter are expected to be intimately related to a diverging ghost 

propagator. In addition, iri^we have reported on a first SU (3) lattice computation of the 
ghost-gluon vertex at zero gluon momentum. We confirm, what has been found already 

for SU(2), namely that the data are quite consistent with a constant vertex functionES. 

For the given report we have extended our investigations restricted so far to pure 
SU(3) Yang-Mills theory to the full QCD case with Nf = 2 clover-improved Wilson 
fermions. For the latter case we have used lattice field configurations produced by the 
QCDSF collaboration and made available via the International Lattice DataGrid (ILDG) b . 



2 Landau gauge fixing, lattice gluon and ghost propagators 

As a first step we have studied the gluon and ghost propagators in the quenched approx- 
imation. SU(3) gauge field configurations U = {U X}fl } thermalized with the standard 

a We thank A. Lokhov for bringing us the arguments iri^to our attention. 

6 We thank Gerrit Schierholz and the QCDSF members for giving us access to their configurations. 



Wilson gauge action have been put into the Landau gauge by iteratively maximizing the 
gauge functional 
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with g x G SU(3). It has numerous local maxima (Gribov copies), each satisfying the 
lattice Landau gauge condition 

(d/A^x) = 9 A^x + ft/2) - %(x - A/2) = (4) 

for the gauge transformed lattice potential 



, • (5) 

traceless 



To explore to what extent this ambiguity has a significant influence on gauge depen- 
dent observables, we have gauge-fixed each thermalized configuration a certain number 
of times (several tens depending on lattice size and coupling) mostly using the over- 
relaxation algorithm, starting from random gauge copies. For each configuration U, we 
have selected the first (f c) and the best (be) gauge copy (that with the largest functional 
value) for subsequent measurements. We have checked that for the be copy the number 
of trials was sufficient for a convergence of the ghost propagator within the statistical 



noise. For details we refer to 
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It turned out that the Gribov ambiguity has no systematic influence on the infrared 
behaviour of the gluon propagator. This holds as long as one restricts to periodic gauge 

transformations on the four-torus. In the gauge orbits have been extended to the full 
gauge symmetry including also non-periodic Z(JV) transformations. As a consequence 
the Gribov problem appears to be enhanced even for the gluon propagator. Here we 
restrict ourselves to the standard case of periodic gauge transformations. In marked con- 
trast to the weak Gribov copy dependence of the gluon propagator the ghost propagator 
at lower momenta depends on the selection of gauge copies. (For a study of the influence 

of Gribov copies on the ghost propagator in the SU(2) case se 

The gluon propagator D^ u (q 2 ) is the Fourier transform of the gluon two-point func- 
tion, i.e. the colour-diagonal expectation value 

DfM = (K{k)A b u {-k)) (6) 

where A®(k) denotes the Fourier transform of A^(x + /t/2) and q is the physical mo- 
mentum 

related to the integer-valued lattice momentum 6 (—L u /2, L^/2] for the linear lat- 

tice extension L^,// = 1, ... ,4. According to Ref. , a subset of possible lattice mo- 
menta k has been chosen for the final analysis of the gluon propagator, although the Fast 



Fourier Transform algorithm provides us with all lattice momenta. In what follows for 

the quenched case we use the lattice spacing a in physical units as determined inPH 

The ghost propagator is derived from the Faddeev-Popov (F-P) operator, the Hessian 
of the gauge functional given in Eq. (|3j. We expect that the properties of the F-P operator 
differ for the different maxima of the functional (Gribov copies). This has consequences 
for the ghost propagator as is shown below. The F-P operator can be written in terms of 
the (gauge-fixed) link variables U Xjfl as 
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2 ReTr 



{T a ,T b }(U x ^ + U x ^) 
T b T a U X J , 



T a T b U x . 



and T a , a = 1,...,8 being the (hermitian) generators of the su(3) Lie algebra 
satisfying Tr [r a T 6 ] = 5 ab /2. The ghost propagator is then determined by inverting 
the F-P operator M 

G ab (q) = ^(e-^^^fM- 1 ]^)^, (9) 

= 5 ab G(q 2 ). 



Following RefJ ' I we have used the conjugate gradient (CG) algorithm to invert M on 
a plane wave ip c with colour and position components ip^(x) = 5 ac exp(27ri k-x). In fact, 
we applied the pre-conditioned CG algorithm (PCG) to solve M"^ b (y) = ip^(x). As 

the pre-conditioning matrix we used the inverse Laplacian operator A -1 with diagonal 
colour substructure. This has significantly reduced the required amount of computing 

time (for details see^j). After solving M(f = tp c the resulting vector <p is projected 

back on ip c such that the average G cc {q) over the colour index c (divided by V) can be 
taken explicitly. Since the F-P operator M is singular if acting on constant modes, only 
k ^ (0, 0, 0, 0) is permitted. Due to high computational requirements to invert the F-P 
operator for each k, separately, the estimator on a single, gauge-fixed configuration is 
evaluated only for a preselected set of momenta k. 



3 Quenched QCD results 

In Fig. ^ we show the ghost and the gluon dressing functions Zg,d(q 2 ) versus q 2 for 
f c copies in the present state of our quenched QCD simulations. Both propagators have 
been renormalized separately for each j3 with the normalization condition Zq n = 1 at 
q = A GeV. 
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Figure 1: The dressing functions for the ghost Zq = q 2 G(q 2 ) (l.h.s.) and gluon propagator Zo 
q 2 D(q 2 ) (r.h.s.) vs. q 2 for quenched QCD, both measured on f c gauge copies. 
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Figure 2: The ratios Z ic /Z hc for the dressing functions Zg (l.h.s.) and Zb (r.h.s.) determined on first (f c) 

and best (be) gauge copies vs. momentum q 2 . 



In Fig. |2] we illustrate the effect of the Gribov copies for periodic gauge transfor- 
mations. We have plotted some fc - to - be ratios of the ghost and gluon dressing 
functions. Obviously there is no influence visible for the gluon propagator within the 
statistical noise. On the contrary, for the ghost propagator the Gribov problem can cause 
0(5%) deviations in the low-momentum region (q < 1 GeV). For better gauge copies 
the ghost dressing function becomes less singular in the infrared. A closer inspection 
of the data for the ghost propagator indicates that the influence of Gribov copies be- 
comes weaker for increasing lattice size. This is in agreement with a recent claim by 
Zwanziger according to which in the infinite volume limit averaging over gauge copies 
in the Gribov region should lead to the same result as averaging over copies restricted to 

the fundamental modular regioiJ^U IiJ^I similar indications have been found for SU (2) 
taking non-periodic Z(2) transformations into account. 

Note in Fig. Q the deviation of the lowest-momentum Zq data point obtained on an 
asymmetric lattice 24 3 x 48. This deviation deserves further study and sheds some critical 



light on investigations on strongly asymmetric lattices^SESl Still it remains difficult to 
extract a possible power behaviour at low q 2 as in Eq. Q. 

In Fig.|3]we present the combined result for the running coupling according to Eq. ^ 
for f c copies in the quenched QCD case. Fits to the 1-loop and 2-loop running coupling 
are also shown. Our data confirm the running coupling monotonously to decrease with 
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Figure 3: L.h.s.: The momentum dependence of the running coupling a s (q 2 ) for quenched QCD measured 
on f c gauge copies. R.h.s.: The inverse ghost-gluon vertex renormalization function Z^ 1 (q 2 ) measured on 
f c gauge copies and normalized to unity at q — 4 GeV. 

decreasing momentum in the range q 2 < 0.3 GeV 2 . The nice coincidence of the data for 
various lattice sizes (except the points refering to the asymmetric lattice) makes it hard to 
see how finite- volume effects can be blamed for the disagreement with the DS continuum 
result Eq. 0. 

One can ask, whether the ghost-gluon vertex renormalization function Z\ (q 2 ) is really 
constant at lower momenta. A recent investigation of this function defined at vanishing 

gluon momentum for the SU(2) case^^ supports that Z\{q 2 ) w 1 at least for momenta 
larger than 1 GeV. We have performed an analogous study for Z\(q 2 ) in the case of 
SU (3) gluodynamics. Our results are presented on the right hand side of Fig.|3] There is a 
slight variation visible in the interval 0.3 GeV 2 < q 2 < 5 GeV 2 . But this weak deviation 
from being constant will not have a dramatic influence on the running coupling. There is 
a data point at the lowest available momentum derived from the asymmetric lattice which 
deviates in the opposite direction. Current simulations on larger lattices will enable us to 
draw conclusions in the near future. 



4 Full QCD results 

Now let us turn to the full QCD case. We have investigated ILDG gauge field configu- 
rations generated in simulations by the QCDSF collaborationist^. jVy = 2 flavours of 
dynamical clover-improved Wilson fermions together with the Wilson plaquette gauge 
action have been used. The (asymmetric) lattice size is 24 3 x 48 possibly demanding 
some caution in the very infrared as discussed already in the quenched case. The values 



Table 1: The parameter values /3, k etc. used for our investigation. 
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5.29 


0.13550 


0.136410(09) 


0.0246 


2.197 


60 


5.29 


0.13590 


0.136410(09) 


0.0138 


2.324 


55 


5.25 


0.13575 


0.136250(07) 


0.0135 


2.183 


60 



of the bare coupling (3 and of the quark mass ma 



1 (\ 1 
ma = - I 



(10) 



taken into account are collected in Tabled The values for k c can be found in™ The last 
column of the Table gives the number of configurations investigated. The lattice spacing 

has been translated into physical units via the Sommer scale r$/a as determined inl^J 
with ro = 0.5 fm. For fixing to Landau gauge we have employed the over-relaxation or 
the Fourier-accelerated gauge-fixing method. Again we concentrate on f c copy results 
postponing the question of the influence of Gribov copies in the full QCD case to a future 
publication. 

In Fig. [4] we present our full QCD results for the ghost and gluon dressing functions 
and, for comparison, selected data points of the quenched case (i.e. for infinite quark 
mass) for f3 = 6.0 and a 32 4 lattice. All dressing functions have been renormalized such 
that Z G D = 1 at q = 4 GeV. 
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Figure 4: The dressing functions of the ghost propagator Zg (l.h.s.) and the gluon propagator Zd (r.h.s.) 
versus q 2 for full QCD with Nf = 2 flavours of clover-improved Wilson fermions (lattice size 24 3 x 
48). Quenched QCD data for (3 = 6.0 and 32 4 are shown for comparison. The dressing functions are 
renormalized to Z — 1 at q = 4 GeV. All the data refer to f c gauge copies. 

The unquenching effect becomes clearly visible for the gluon propagator, whereas the 
ghost propagator does not show any quark mass dependence. The non-perturbative peak 



of the gluon propagator at q ~ 1 GeV becomes softer as the quark mass is decreasing. 
This has been observed also in recent lattice computations of the gluon propagator using 

dynamical AsqTad improved staggered quarksl^as well as fro m unqenchi ng studies for 

the ghost and gluon propagators within the DS equation approach P^^^ ! We refer also 

to recent lattice studies with dynamical Kogut-Susskind and Wilson fermions^OJ^JJ 

In Fig. |3] we show the corresponding MOM scheme running coupling a s {q 2 ) and the 
inverse of the renormalization function Z\(q 2 ) of the ghost-gluon vertex as functions 
of the physical momentum. Again we compare results for full QCD with those for the 
quenched case at (3 = 6.0 on a 32 4 lattice. The normalization factors of the running 
couplings were determined by fitting the data to the (massless) 2-loop coupling formulae 
for a s (q 2 ) with Nf = 2 and Nf = 0, respectively, both in the range q 2 > 30 GeV 2 . The 
vertex renormalization function was normalized to Z\ = 1 at q = 4 GeV for each data 
set (P, ma) separately. 

For the running coupling we observe a quite clear unquenching effect extending from 
the perturbative range down to the infrared region. Again we see the coupling to turn 
down in the infrared limit. However, the data points at q 2 ~ O.lGeV 2 have to be checked 
carefully on symmetric and larger lattices. The same holds for Z\ . 
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Figure 5: The running coupling (l.h.s.) and the inverse renormalization function Z^ 1 of the ghost-gluon 
vertex (r.h.s.) are shown as functions of the physical momentum for full QCD. For comparison, the quenched 
QCD running coupling data for (3 — 6.0, 32 4 are shown with a spline fit line to guide the eye. The 2-loop 
fit to a s (q 2 ) for Nf = 2 is drawn with the solid line. The vertex renormalization function is normalized to 

Zi = 1 at q = 4 GeV. 



5 Conclusions 

We have studied the low-momentum region of QCD within the Landau gauge using 
Monte Carlo simulations with the Wilson plaquette action and various lattice sizes from 
16 4 to 32 4 for the quenched case in the range of bare couplings = 5.8, . . . , 6.2. In 
this way we have reached momenta down to q 2 ~ O.lGeV 2 . Moreover, for studying 



unquenching effects we have carried out an analogous investigation in full QCD for the 
case of two improved Wilson flavours for a couple of quark mass values but only one 
lattice size 24 3 x 48. Our data presented here refer to first gauge copies as obtained from 
over-relaxation or Fourier-accelerated gauge fixing. In the quenched QCD case we have 
seen clear Gribov copy effects in the infrared region for the ghost propagator. The ghost 
propagator becomes less singular within a 0(5%) deviation, when better gauge copies 
are taken. We have found indications that the Gribov effect also weakens as the volume 
is increasing. Towards q — > in the infrared momentum region, the gluon dressing func- 
tions in quenched as well as full QCD were shown to decrease, while the ghost dressing 
functions turned out to rise. However, the connected power laws predicted by the infinite- 
volume DS approach could not be confirmed on the basis of our data. Correspondingly, 
the behaviour of the running coupling a s (q 2 ) in a suitable momentum subtraction scheme 
(based on the ghost-gluon vertex) does not approach the expected finite limit. Instead, 
this coupling has been found to decrease for lower momenta after passing a turnover at 
q 2 w 0.4 GeV 2 for quenched (N f = 0) as well as for full QCD (N f = 2). Unquenching 
effects have been clearly recovered for the gluon propagator, whereas the ghost propa- 
gator was almost unchanged. This is in one-to-one correspondence with what has been 
recently found in the Dyson-Schwinger equation approach. However, the puzzle of the 
existence of a non-trivial infrared fixed point in the infinite volume limit remains un- 
solved. 
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